Draft version October 8, 2009 

Preprint typeset using I4TgX style emulateapj v. 08/22/09 



HALO RETENTION AND EVOLUTION OF COALESCING COMPACT BINARIES IN COSMOLOGICAL SIMULATIONS 
OF STRUCTURE FORMATION: IMPLICATIONS FOR SHORT GAMMA-RAY BURSTS 

Marcel Zemp 1 ' 2 , Enrico Ramirez-Ruiz 2 and Jurg Diemand 2,3 

Draft version October 8, 2009 

ABSTRACT 

Merging compact binaries are the one source of gravitational radiation so far identified. Because short-period 
systems which will merge in less than a Hubble time have already been observed as binary pulsars, they are 
important both as gravitational wave sources for observatories such as LIGO but also as progenitors for short 
gamma-ray bursts (SGRBs). The fact that these systems must have large systemic velocities implies that by the 
time they merge, they will be far from their formation site. The locations of merging sites depend sensitively 
on the gravitational potential of the galaxy host, which until now has been assumed to be static. Here we 
refine such calculations to incorporate the temporal evolution of the host's gravitational potential as well as 
that of its nearby neighbors using cosmological simulations of structure formation. This results in merger site 
distributions that are more diffusively distributed with respect to their putative hosts, with locations extending 
out to distances of a few Mpc for lighter halos. The degree of mixing between neighboring compact binary 
populations computed in this way is severely enhanced in environments with a high number density of galaxies. 
We find that SGRB redshift estimates based solely on the nearest galaxy in projection can be very inaccurate, 
if progenitor systems inhere large systematic kicks at birth. 

Subject headings: gamma rays: bursts — stars: formation — cosmology: observations — galaxies: formation 
— methods: N-body simulations 



1. INTRODUCTION 

The association of short gamma-ray bursts with both star- 
forming gala xies and with ellipticals dom inated by old stellar 
populations (|Berger et al. 2005 ; Bloom et al. 200 6~l|Fox et al.| 
[2^03|[Gehrels et al.|2005||prochaska et al.|2006||Berger|20091 l 
suggested an analogy to type la supernovae, as it indicated a 
class of progenitors with a wide distribution of delay times 
between formation and explosion. Similarly, just as core- 
collapse supernovae are discovered almos t exclusively in star- 
forming galaxies, so too are long GRBs (Woosley & Bloom 
2006). Indeed, a detailed census of the types of host galaxies, 
burst locations and redshifts could help decide between the 
various SGRB p rogenitor alternatives (e. g. Zheng & Ramirez- 
|Ruiz|2007|rGuetta & Piran|2005||BToom & Prochaska 2006): 
If the progenitor lifetime is long and the systemic kick is 
small, then the bursts should correspond spatially to the old- 
est populations in a given host galaxy. For early-type galax- 
ies, the distribution would most likely follow the light of the 
host. A neutron star (NS) binary could take billions of years to 
spiral together, and could by then, if given a substantial kick 



velocity on formation, merge far from its birth site (|Fryer et 



al. 1999; B loom et al.|1999[|Belczynski et al.|2002||Bloom et 
al.|2002|) . The burst offsets would then most likely be larger 
for smaller mass hosts. 

Double neutron star binaries, such as the famous 
PSR1913+16, will eventually coalesce, when gravi tational ra- 
diation drives them together ( Kalo gera et al. [2007) . Each su- 
pernova is thou ght to impart a sub stantial kick to the result- 
ing NS (Hansen & Phinney 1997). For systems that survive 
both supernovae explosions the center of mass of the rem- 
nant binary itself will receive a velocity boost on the order 
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of a few hundred kilometers per seco nd ( |Brandt & Podsiad- 
|lowski|1995"t [Fryer & Kalogera|1997] >. As a result, NS bina- 
ries will be ejected from their birth sites. The exact distribu- 
tion of merger sites depends sensitively on the gravitational 
potential of the host, which until now has been assumed to be 
static. The potential of a realistic host galaxy is, however, not 
static. In fact, the gravitational potential of the host as well as 
that of its nearby neighbors is expected to evolve dramatically 
from compact binary production until coalescence. In order to 
incorporate these effects self-consistently, in this Letter, we 
study the orbital evolution of compact binary systems using 
cosmological simulations of structure formation. Our results 
provide new insights into what happens when compact binary 
stars are ejected from their birth halos as a result of veloc- 
ity kicks, and what progenitor clues a distant observer might 
uncover from the distribution of SGRB sites in and around 
galaxies. 

2. COMPACT BINARIES IN COSMOLOGICAL SIMULATIONS 

The focus of this work is to understand the retention and 
evolution of compact binaries in an evolving cosmological 
simulation. To this end, we have performed a dark matter 
only cosmological structure formation simulation. A 80 co- 
moving Mpc periodic box with a single mass resolution of 
m p = 1.07 x 10 9 Mq corresponding to 256 3 particles is initial- 
ized at a starting redshift of z = 22.4 (161 Myr). The parti- 
cles have a softening length of 16 kpc a nd we have used the 
WMAP 3-year cosmological parameters (Sp ergel et al.|2007) 
with fi M ,o = 0.238, ft A ,o = 0.762 and H = 73 km s _1 Mpc" 1 , 
erg = 0.74 and n s = 0.951. The time evolution was done 
with the parallel tree-code PKDGRAV2 ( |Stadel|200l) . PKD- 
GRAV2 uses a fast multipole expansion technique in order to 
calculate the forces with a hexadecapole precision and a time- 
stepping scheme that is based on the true dynamical t ime of 
the particle s with an accuracy parameter of i]£, = 0.03 (Zemp 
|et al.|2007) . 
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First, we evolve the initial conditions until redshift z = 1 .60 
(4.24 Gyr). At this time we find all the halos in our sim- 
ulation with a friends-of-friends method ( |Davis et al.||1985] l 
and select all those with a minimum of 200 particles of which 
there are 2461. This criterion corresponds to a halo mass of 
2.15 x 10 11 Mq. Second, we populate each of these selected 
halos with 2000 massless tracer particles which are placed 
at the centre of each halo. Each tracer particle is meant to 
represent a compact binary system, which, on a verage, forms 
aroun d the peak of the star formation epoch (Madau et al. 
1996). The velocity distribution of the tracer particles is as- 
sumed to be isotropic and to have a Maxwell-Boltzmann dis- 
tribution with a mean speed of v = 360 km s -1 and a disper- 
sion of a = 150 km s -1 . This is consistent with the magni- 
tude of the natal kicks required to explain the observed pa- 
rameters of binary neutron star systems - only when kicks 
have magnitudes exceeding 200 km s -1 can the progenitor or- 
bits be sufficiently wide to accommodate evolved helium stars 
and still produc e the small separations measured in these sys - 
tems ( |Brandt & Podsiadlowski|1995||Fryer & Kalogera|1997] l. 
Third, we weight each tracer particle by w, = m,/m max where 
nti is the mass of the halo where the tracer particle with in- 
dex i initially starts and m max is the mass of the most massive 
halo at redshift z = 1 .60. This is done in order to account for 
the contribution of a given halo to the total number of com- 
pact binaries produced at z = 1.60 under the assumption that it 
scales with the mass content in each halo. Finally, we evolve 
the cosmological cube together with the tracer particle popu- 
lations until redshift z = (13.8 Gyr). 

To explore the ability of individual halos in retaining their 
own birth compact binary population, one requires to accu- 
rately follow not only the location of the tracer particles but 
also the fate of the individual halos as they evolve and expe- 
rienced a substantial metamorphosis. To this end, we have 
marked the particles with the highest phase space density in 
each of the 2461 halos that have been populated with tracer 
particles. By following these marked particles, we can then 
accurately track the location of the halos over time. These 
high phase space density particles are optimal for tracking the 
halo positions since they stay at the centre of the halo and are 
only minimally affected by tidal effects. Even when a halo be- 
comes tidally disrupted, the centre of the debris of these high 
phase space density particles still provides a good represen- 
tation of where the halo would be if it would not have been 
disrupted. 

3. THE COSMOLOGICAL COMPACT BINARY MERGING SITES 

We now turn to an examination of the cosmological distri- 
bution of the compact binary systems, which is done here by 
following the massless tracer particle populations over time. 
Herefore, we picked out by eye three halos, each in a differ- 
ent environment (field, group and cluster). Figure [T] shows 
the probability density function of finding tracer particles as a 
function of distance from a given halo at three different snap- 
shots in time. The red lines show the radial probability density 
function for all tracer particles originating in the particular 
halo while the blue lines show the same function for tracers 
belonging to all other halos. 

Figure [TJ is self-explanatory. In the field environment, the 
distribution is dominated out to a distance of a few Mpc by its 
own tracer particles at all cosmological times. The tracer par- 
ticle population injected at z = 1.60 (4.24 Gyr) is observed to 
be rather extended ~ 10 Mpc. This is because this field halo, 



whose mass at z = 1.60 was 2.15 x 1O H M , was unable to 
effectively retain most of its own tracer particles. This is also 
the case for the host halo belonging to the group environment, 
whose central mass at z = 1.60 was 2.27 x 10 n M Q . However, 
in this more crowded environment, the close proximity of the 
neighboring halos allows foreign tracer particles to pollute the 
central regions of the host halo. In the cluster environment, 
only very few tracer particles are able t o escape the deep halo 
potential well where they were born (Niino & Totani 2008), 
whose mass at z = 1.60 was 5.95 x 1O^M . The few unbound 
particles are still effectively retained by the cluster's potential. 
As a result of the high merger activity in such cluster environ- 
ments, the mixing of the various tracer populations increases 
dramatically with time. At z = 0, for example, the probabil- 
ity of finding a foreign tracer is equal or higher at all radial 
distances than finding one originating from the massive cen- 
tral halo, whose mass at z = is now 7.64 x 10 14 M Q . This 
clearly indicates, that in a high density environment with a lot 
of merger activity there is a high degree of tracer mixing and, 
as a result, the closest galaxy at the time of merging is likely 
not to be the one where the compact binary system originated 
from. 

The ability of a halo to retain its birth population of tracer 
particles is further illustrated in Figure [2] where we plot the 
density field at z = (left column) of the three halo environ- 
ments shown in Figure [T] together with the distribution of the 
particle tracers that originate from the selected halo (right col- 
umn). We see that in the cluster environment, the tracer par- 
ticle population stays compact. In this crowded environment, 
the many accreted subhaloes contribute significantly to the 
pollution of the tracer particle population in the central clus- 
ter region as shown in Figure [T] In the small mass field halo, 
the tracer particles spread over many Mpc in distance. Given 
that there are not many other halos in the immediate neigh- 
bourhood, the tracer particle pollution is negligible. This is, 
however, not the case for the selected halo in the group en- 
vironment which despite having a mass similar to that of the 
field halo at z = 1 .60, experiences a substantial degree of tracer 
particle mixing from the other group members. 

Now the question arises over whether the compact binary 
was born in the galaxy it is closest to at the time of merger? 
If SGRB are triggered by compact binary formation, then the 
merging site can be deduced by the afterglow location. We 
can answer this by following the population of compact bi- 
naries as the cosmological simulation evolves. At any given 
time we can calculate the weighted fraction of tracer particles 
that are still closest to the halo they started at redshift z = 1 .60 

given by / = W k/T,f m Wu where N tol = 2461 x 2000 = 
4922000 is the total number of tracer particles and the k sum- 
mation runs over all N co tracers that are still closest to their 
original halo. In Figure [3] we plot this fraction as a function 
of the environment mass containing the tracer particle popu- 
lation. 

The environment of a tracer particle is defined by the most 
massive halo we find at a given time that contains the tracer 
particle. We define that a tracer particle is contained within a 
halo if they are within a sphere of radius d = 2s from its centre 

and where s is given by s = \ (j^Li **^-**,minJ and x ktIaax 

respectively -x>t. m in are the maximum/minimum coordinate in 
the kth dimension of any particle in the halo. The fraction / 
is not very sensitive to the detailed definition of environment 
used here since for example the definition d = s results in a 
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qualitatively similar plot. 

Figure [^highlights the importance of environment on com- 
pact binary retention. It shows that in environments < 
10 12 M Q around 95% of all weighted tracers at all cosmo- 
logical times are still closest to their halo where they origi- 
nate from. For more massive environments, this ceases to be 
true. For example, in environments > 10 14 M & , the fraction 
of binaries found closest to their birth halo severely decreases 
with time, reaching 45% at redshift z = 0. This implies that 
in more than half the cases the merger event would not take 
place closest to its birth galaxy site if occurring within a clus- 
ter environment. 

4. DISCUSSION 

Before the detailed modeling of light curves was used to 
constrain the nature of supemovae progenitors, the location of 
supernovae in and around galaxies provided important clues 
to the nature of the progenitors. Si milarly, in th e absence of 
supernova-like features (e.g. Hjorth et al. 2005a), detailed ob- 
servations of the astrophysics of individual host galaxies may 
thus be essential before stringent constra ints on the identity 



of SGRB progenitors can b e placed 4 (e.g. Zheng & Ramirez- 
|Ruiz|2007[|Fox et al |2005] l. 

Even with a handful of SGRBs detected to date, it has be- 
come apparent that short and long events are no t drawn from 
the same parent stellar population (Nakar 2007). In contrast 
to long GRBs, the galaxies associated with SGRBs exhibit a 
wide rang e of star-fo rmation rates, morphologies and metal- 
licities (Berger 2009). They are also often found in older and 
lower-redshift galaxies and, in a few cases, with large (> 10 
kpc) projected offsets from the centers of their putative host 



galaxies (Bloom et al. 2006, Prochaska et al. 2006; Berger et 
2005j [Gorosabel et al.|2006t |Levan et al.|2006|[Bloom et 



al. 



al. 2001) . SGRBs are, however, not universally at large offsets 
and are not always associated with early-type galaxies (|Fox et 
al.|2005[|ViTiasenor et al.|2005l|Hjorfh et al.|2005b||Covino et 
al. | 2006 |Sodeiter g"et^rpn06| l. 



The discovery of early-type galaxy hosts suggest a progeni- 
tor lifetime distribution extending well beyond a Gyr. A large 
progenitor lifetime would help explain t he apparent high inci- 
dence of galaxy cluster membership (|Bloom e t al. 2006] [Ped^] 
|ersen et aL"1|2005| |Berger et al.||2007a| l. On the other hand, 
shorter lifetimes are required to explain the population of 
SGRB at moder ately high redshift (Berger et al.|2007b||Gra~ 



ham et al. 2009 ). The observed projected distances from what 



has been argued are the plausible hosts, if true, also ho lds im- 
portant ramifications for the sort of viable progenitors (Bloom 
|& Prochaska |[200 6). The large offsets seen from early- type 
hosts would seem to be at odds with progenitor syste ms with 
small systematic kicks (such as in globular cluste rs Grindlay 
[etaLl[2006l [Lee & Ramirez-Ruiz|2007l |Parsons et al.|2009^ , 
although with such large physical offsets the possibility re- 
mains that the association with the putative host is coinciden- 
tal. On the other hand, ba sed on the small offsets from some 
lo w-mass galaxy hos ts (Prochaska et al. 2006; Soderberg et 
aLl |2006l |Berger|[2009) |Bloom & Prochaska ||2006| , SGRB 
progenitors cannot all have large systematic kicks at birth and 
inhere large delay times from formation. 

Making more quantitative statements about the nature of 
the progenitor systems is not only hampered by small number 
statistics but also from the lack of robust predictions of the 
distribution of merger sites. These distributions depend sen- 



sitively on the gravitational potential of the host, which until 
now has been assumed to be static, and the compact binary 
formation properties, especially the systematic kick velocity. 

Here we have refined such calculations to incorporate the 
temporal evolution of the host's gravitational potential as well 
as that of its nearby neighbors self-consistently using cosmo- 
logical simulations of structure formation. This results in di- 
verse predictions of offsets and compact binary demographics 
even in the simplest case of a kick velocity distribution (here 
assumed to be in excess of 200 km s" 1 in order to expl ain the 
observed pa rameters of double neutron star systems; |Fryer & 
Kalogera 1997 ) whose properties do not vary with the initial 
binary separation. 

Two important predictions stand out. First, the merger site 
distributions computed in this way are more diffusively lo- 
cated with respect to their putative hosts. In a field environ- 
ment, for example, the distribution of merging sites can ex- 
tend out to a distance of a few Mpc. This is more severe for 
those host galaxy halos that were unable to effectively retain 
most of its own compact binary population at birth. Second, 
the degree of mixing between neighboring compact binary 
populations depends on galactic environment. In a cluster, 
for example, the mixing of the various compact binary pop- 
ulations is severe as a result of the high merger activity and 
increases dramatically with time. At z = 0, in a cluster envi- 
ronment, the probability of finding a foreign coalescing com- 
pact binary system is equal or higher at all radial distances 
than finding one originating from the massive central halo. As 
a result, the closest galaxy at the time of binary coalescence 
and possibly SGRB occurrence may not to be the one where 
the compact binary system originated from. 

Of course, our basic model is rather simple since we assume 
a single epoch of star formation and a simple star formation 
recipe. Also different distributions of kick velocities should 
be considered. We do not expect that our qualitative results 
will change dramatically as our first, more general, results in- 
dicate. 

It is evident form the discussion above that assuming a 
large (already evolved) host galaxy at the time of compact 
binary formation thus severely overestimates the binary re- 
tention fraction and the concentration of their merging site 
distribution. This implies that SGRB redshift estimates based 
solely on the nearest galaxy in projection can be very inac- 
curate, if progenitor systems inhere large systematic kicks at 
birth. Interpretations on the nature of the SGRB progenitor 
using the stellar and mass properties of the nearest galaxy in 
projection as established by the afterglow location must there- 
fore be regarded with suspicion. Finally, it should be noted 
that a direct comparison with model predictions is s till im- 
peded by the possibility o f an ambient density bias (Bloom 
|et al.|2006[|Lee et al.|2005] l where SGRBs are more likely to 
be found in denser gas regions and , as a result, we could be 
missing a population of bursts with large systematic kicks. 



We thank J. Bloom, S. Faber, C. Fryer, V.Kalogera, R. 
O'Shaughnessy and X. Prochaska for useful discussions. 
M.Z. is supported by NSF grant AST-0708087. We fur- 
ther acknowledge support from Swift: NASA NNX08AN88G 
(E.R.), NSF: 0521566, and the David and Lucile Packard 
Foundation (E.R.). 



4 



Zemp, Ramirez-Ruiz & Diemand 



4 The interaction of burst ejecta with a stell ar binary companion {Mac-| 
Fadyen et al. 2005 1 or with its emitted radiation I Ramirez-Ruiz 2004 1 could 

REFE1 

Belczynski, K., Bulik, T., & Kalogera, V. 2002, ApJ, 571, L147 
Berger, E., et al. 2005, Nature, 438, 988 

Berger, E., Shin, M.-S., Mulchaey, J. S., & Jeltema, T. E. 2007, ApJ, 660, 
496 

Berger, E., et al. 2007, ApJ, 664, 1000 
Berger, E. 2009, ApJ, 690, 231 

Bloom J. S., Sigurdsson S., Pols O. R., 1999, MNRAS, 305, 763 
Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 1111 
Bloom, J. S., et al. 2006, ApJ, 638, 354 

Bloom, J. S., & Prochaska, J. X. 2006, in AIPCS, ed. S. S. Holt, N. Gehrels, 

& J. A. Nousek, 836, 473 
Bloom, J. S., et al. 2007, ApJ, 654, 878 
Brandt, N., & Podsiadlowski, P. 1995, MNRAS, 274, 461 
Covino, S., et al. 2006, A&A, 447, L5 
Davis, M., et al. 1985, ApJ, 292, 371 
Fox, D. B., et al. 2005, Nature, 437, 845 
Fryer, C, & Kalogera, V. 1997, ApJ, 489, 244 

Fryer, C. L., Woosley, S. E., & Hartmann, D. H. 1999, ApJ, 526, 152 
Gehrels, N., et al. 2005, Nature, 437, 851 
Gorosabel, J. et al. 2006, A&A, 450, 87 

Grindlay, J., Portegies Zwart, S., & McMillan, S. 2006, Nature Physics, 2, 
116 

Graham, J. E, et al. 2009, ApJ, 698, 1620 

Guetta, D., & Piran, T. 2005, A&A, 435, 421 

Hansen, B. M. S., & Phinney, E. S. 1997, MNRAS, 291, 569 

Hjorth, J., et al. 2005a, ApJ, 630, LI 17 



also help shedding light on the identity of the progenitor system. 



Hjorth, J., et al. 2005b, Nature, 437, 859 

Kalogera, V., Belczynski, K, Kim, C, O' Shaughnessy, R., & Willems, B. 

2007, Phys. Rep., 442, 75 
Lee, W. H, Ramirez-Ruiz, E., & Granot, J. 2005, ApJ, 630, L165 
Lee, W. H, & Ramirez-Ruiz, E. 2007, New Journal of Physics, 9, 17 
Levan, A. J., et al. 2006, MNRAS, 368, LI 

MacFadyen, A. I., Ramirez-Ruiz, E., & Zhang, W. 2005, astro-ph/05 10192 
Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C., 

& Fruchter, A. 1996, MNRAS, 283, 1388 
Nakar, E. 2007, Phys. Rep., 442, 166 
Niino, Y., & Totani, T. 2008, ApJ, 677, L23 

O' Shaughnessy, R., Belczynski, K, & Kalogera, V. 2008, ApJ, 675, 566 

Parsons, R. K, Ramirez-Ruiz, E., & Lee, W. H. 2009, arXiv:0904.1768 

Pedersen, K, et al. 2005, ApJ, 634, L17 

Prochaska, J. X., et al. 2006, ApJ, 642, 989 

Ramirez-Ruiz, E. 2004, MNRAS, 349, L38 

Soderberg, A. M., et al. 2006, ApJ, 650, 261 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Stadel, J. G. 2001, Ph.D. Thesis, University of Washington 

Villasenor, J. S., et al. 2005, Nature, 437, 855 

Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 

Zemp, M., Stadel, J., Moore, B., & Carollo, C. M. 2007, MNRAS, 376, 273 

Zheng, Z., & Ramirez-Ruiz, E. 2007, ApJ, 665, 1220 



Halo Retention and Evolution of Coalescing Compact Binaries 5 



field environment group environment cluster environment 

I i i 1 1 mij i i i iiiiij mi iiiiij i i ■ mi iiiiij i 1 1 ■ 1 1 1 ^ i i 1 1 mij i i m i in inij i i i niiij i i i iiiiij ill 




r [kpc] 

FIG. 1 . — Probability density function P(r) of the tracer particles as a function of distance from a particular halo for three different environments and at three 
different snapshots in time. The dash-dotted red lines show the radial probability density function for all tracer particles formed in the particular halo while the 
solid blue lines show the same function for tracers belonging to all other halos. Note how the bound and unbound tracers of local origin (dash-dot/red) give rise 
to bimodal distributions in the field and group environments, while all of them remain bound in cluster environment. 
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FIG. 2. — The tracer particle populations as a function of environment in the local universe. The dark matter density field at z = (top row) of the three halo 
environments shown in FigureMlare plotted together with the distribution of the particle tracers that originate from the selected central halo (bottom row). The 
side length of each panel is 10 Mpc. 
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FIG. 3. — The weighted fraction / of tracer particles over time that are still closest to the halo they started at redshift z = 1.60 as a function of the mass of the 
most massive nearby halo (see text for details). 



